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Using data from 2.9 fb -1 of integrated luminosity collected with the CDF II detector at the 
Tevatron, we search for resonances decaying into a pair of on-shell gauge bosons, WW or W Z, 
where one W decays into an electron and a neutrino, and the other boson decays into two jets. 
We observed no statistically significant excess above the expected standard model background, and 
we set cross section limits at 95% confidence level on G* (Randall-Sundrum graviton), Z' , and W' 
bosons. By comparing these limits to theoretical cross sections, mass exclusion regions for the three 
particles are derived. The mass exclusion regions for Z' and W' are further evaluated as a function 
of their gauge coupling strength. 
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Some models of new physics beyond the standard 
model predict particles that decay into pairs of on-shell 
bosons, for example Z' , or the Randall-Sundrum 

graviton G* 0). Searches for these particles in different 
decay channels have been reported elsewhere [3rtZ[ . Most 
of them used final states consisting of only leptons or 
photons. In this Letter we search for these particles in 
the form of diboson resonances where one boson is a If 
decaying into an electron and a neutrino, and the other 
is a W or Z which decays into two jets. This search 
has the advantage of detecting two types of diboson res- 
onances, WW and WZ, with the same final-state topol- 
ogy. The hadronic decay mode of the W or Z to two jets 
has a higher branching fraction compared to the leptonic 
mode; however, the background from jets also increases. 
Thus we implement a selection based on transverse en- 
ergy (Et) [8j of the detected objects in the final state to 
reduce standard model backgrounds and enhance sensi- 
tivity. 

The diboson decay modes of Z' and W' directly probe 
the gauge coupling strength between the new and the 
standard model gauge bosons. The coupling strength 
strongly influences the decay branching ratios and the 
natural widths of the new gauge bosons. In an ex- 
tended gauge model theory [lj the standard model cou- 
pling strength, g ■ cos 6 W , is replaced by £ • g ■ cos W , where 
£ = C ■ (Mw/My) 2 ; C is a parameter that sets the cou- 
pling strength; and My is the mass of the new gauge 
boson, Z 1 or W'. We set cross section limits on Z' and 
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W' as a function of mass and of £. Our results extend 
the sensitivity beyond the CDF Run I W' results Q with 
almost 30 times the integrated luminosity, and, for the 
first time, set Z' limits as a function of mass and gauge 
coupling strength. For G* , the coupling constant k/Mpi 
dictates the branching ratio and natural width Q , where 
k and M pi are respectively the curvature of the extra 
dimension and the reduced Planck mass scale. This is 
also the first search for the 67* in the WW decay mode. 

This analysis is based on data corresponding to an 
integrated luminosity of 2.9 fb _1 collected using the 
CDF II detector between March 2002 and February 2008. 
The detector is approximately forward-backward and az- 
imuthally symmetric. The detector elements relevant to 
this analysis are the tracking system and the calorime- 
ters. The tracking system consists of an eight-layer sili- 
con tracker (loj surrounded by a 96-layer open-cell drift 
chamber (COT) [ll|. The fiducial coverage of the COT 
is 1 77 1 < 1.0 and the silicon detector extends the cov- 
erage to 1 77 1 < 2.0. The integrated tracking system is 
contained within a superconducting solenoid, providing 
a 1.4 T magnetic field. Surrounding the tracking sys- 
tem are the electromagnetic (EM) and hadronic calorime- 
ters [HI, divided into "central" (\r]\ < 1.1) and "plug" 
(1.1 < 1 77 1 < 3.6) regions. The calorimeters arc made of 
lead (EM) and iron (hadronic) absorbers sandwiched be- 
tween plastic scintillators that provide measurements of 
shower energies. At approximately the shower maximum, 
the EM calorimeters contain fine-grained detectors [l3| 
for measuring shower positions and profiles. 

As we are looking for events with an electron, a neu- 
trino, and two jets, we start with data that were collected 
with an online selection requirement of a central elec- 
tron with 1 77 1 < 1 and Et > 18 GeV. From this dataset 
we select events that have an isolated electron [14| with 
Et > 30 GeV, a neutrino identified by the require- 
ment that the missing Et(Et) > 30 GeV, two or three 
jets with \r}\ < 2.5 and Et > 30 GeV, and an overall 
Ht > 150 GeV, where Ht is the scalar sum of the elec- 
tron Et, the Et, and the Et of all jets 15]. 

To form a WW or WZ hypothesis for the selected 
events, the electron and Et are first combined to form a 
W candidate. Because the longitudinal component of the 
neutrino momentum (E%) is not available, the invariant 
mass of the electron and Et is artificially set to the W 
mass. With this assumption, the conservation of energy 
and momentum results in a quadratic equation for E v z . 
If the discriminant of the quadratic equation is negative, 
the combination is discarded. If it is positive, there are 
two solutions and both are kept. In addition, two jets are 
combined to form a second W candidate or a Z candi- 
date. In the case of a W candidate, we require the two-jet 
invariant mass (Mjj) to fall between 65 and 95 GeV/c 2 , 
corresponding to ± 1.5c of the expected reconstructed W 
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resolution. In the case of a Z candidate, this window is 
between 75 and 105 GeV/c 2 . For a three-jet event, there 
are three two-jet invariant mass combinations. In this 
case only the pair with the invariant mass closest to ei- 
ther the W or the Z mass is kept in order to reduce the 
combinatorial background. The reconstructed W or Z 
candidates are then combined to form the final WW or 
W Z invariant mass. 

Twelve standard model processes are considered as 
background for this analysis: W(—t e ± ^)+jets, QCD 
jets, it, WW, Z{-+ e+e-)+jets, W(—> r^+jets, single 
top, WZ, W'y, Z — y t + t~ , 77, and ZZ. The dominating 
background is W+jets whose contribution is estimated 
by Monte Carlo simulation using the ALPGEN jT^| event 
generator, interfaced to pythia [17] for parton shower- 
ing and followed by the geant 3 [l8| based CDF II de- 
tector simulation. With the exception of the QCD jet 
background, the rest of the background processes are all 
estimated by Monte Carlo simulation using the pythia 
event generator. The cross sections used for the simu- 
lated background processes are obtained from NLO cal- 
culations. 

The QCD jet background comes from events with three 
or more jets where one of the jets is misidentified as an 
electron. With this misidentified electron, the event may 
pass through subsequent event selection criteria and the 
reconstruction processes. The contribution of the QCD 
jet background is estimated using a dataset that has an 
online selection requirement of one jet with Et > 20 
GeV. We first exclude events that have any identified 
electrons, then each jet in the central region is treated 
as an electron with a weight corresponding to the prob- 
ability that a jet is misidentified as an electron. This 
probability is a function of jet Et and varies from 10~ 4 
at 30 GeV to 10~ 3 above 100 GeV @. The misidentified 
electron is combined with the $t and then with two jets 
to form WW or WZ candidates as described earlier. The 
resulting QCD jet background is normalized to the data 
by matching the $t spectrum between data and expected 
background at their peaks around 10 GeV, where little 
signal is expected and the QCD jet background domi- 
nates. This normalization factor is used for the QCD 
jet contribution throughout the analysis. Figure [T] shows 
the resulting J£t spectrum for events with two jets that 
would have passed the event selection criteria except for 
the E T > 30 GeV cut. 

The systematic uncertainties taken into account in 
the background calculations are the following, listed by 
decrea sing significance: jet energy scale (JES) uncer- 
tainty [l9(, theoretical cross section uncertainty [201 ], lu- 
minosity uncertainty [2l| . and jet misidentification rate 
uncertainty. The dominating systematic uncertainty is 
the JES uncertainty which amounts to ~ 13% of the es- 
timated background. The cross section and luminosity 
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FIG. 1: Pt spectrum from events with two jets. The vertical 
line marks the fir > 30 GeV cut. The QCD component is 
scaled such that data and expected background match at the 
peak area where no signal is expected. "Others" background 
includes: Z(— > e + e~)+jets, single top, WZ, Wj, Z —y t + t~ , 
77, and ZZ. 

uncertainties are ~ 6% each. 

Signal detection efficiencies arc also determined from 
simulated events using the pythia event generator. For 
a set of selected mass values ranging from 165 GeV/c 2 to 
1000 GeV/c 2 , the three types of particles are simulated: 
G* with k I'M pi = 0.1, Z' and W' with pythia default 
settings corresponding to the extended gauge model with 
a suppression factor £ = (Mw/My) 2 , i.e., C = 1. The re- 
constructed signals are Gaussian in shape, and the mass 
resolution is linearly proportional to the generated mass 
values, varying from 20 GeV/c 2 at 200 GeV/c 2 mass to 
80 GeV/c 2 at 1000 GeV/c 2 mass. For calculating the 
efficiencies we choose an acceptance mass window cor- 
responding to ±1.5 times the reconstructed signal reso- 
lution. This choice gives a good signal to background 
ratio. The same acceptance mass windows are also used 
to obtain the number of background events. 

The systematic uncertainties taken into account for 
the signal acceptance, defined as the product of signal 
detection efficiency and integrated luminosity, in order 
of decreasing significance, are: jet energy scale (JES) 
uncertainty, luminosity uncertainty, initial state radia- 
tion (ISR) uncertainty, final state radiation (FSR) un- 
certainty, and parton distribution function (PDF) uncer- 
tainty. Similarly to the background uncertainties, the 
JES uncertainty dominates the systematic uncertainties 
and varies from 12% at 170 GeV/c 2 mass to 6% at 700 
GeV/c 2 mass for G* , 13% (170 GeV/c 2 ) to 6% (1000 
GeV/c 2 ) for Z' , and 9% (190 GeV/c 2 ) to 6% (1000 
GeV/c 2 ) for W'. ISR, FSR and PDF uncertainties are 
of the order of 1 — 3% each and decrease with increasing 
diboson mass. 

In order to improve sensitivity at higher mass, addi- 
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tional sets of higher Et cuts for the constituent particles 
(observed in the detector as electron, $t from neutrino, 
and jets) are tried. Two series of the Et cut sets are 
implemented. The hrst series requires a higher Et on 
all four participating particles ranging from 40 GeV to 
80 GeV in steps of 10 GeV. The second series requires a 
higher Et on only one daughter particle from each of the 
decaying bosons, i.e., a higher Et for either the electron 
or the neutrino, and the same higher Et for one of the 
two jets. The Et values in this series range from 40 GeV 
to 120 GeV in steps of 10 GeV. For each set of Et cuts 
the systematic uncertainties for the backgrounds and the 
acceptances are re-evaluated, but are found to be not 
very sensitive to the variations. 

To find the optimal set of Et cuts at each selected 
mass point, the expected cross section limits, which are 
based only on the background and the signal acceptance, 
are calculated for each set of cuts. We found that the 
first series of Et cuts gives the best expected limits for 
Z' and W, while the second series is best for G* . The 
optimal Et cuts for each particle type are then selected 
from their own optimal series. The sets that give the 
best expected limits are chosen without reference to their 
impact on the data sample. Although the background 
processes respond differently to the two series of Et cuts, 
the best expected limits obtained from each series are 
very similar. Generally, as the mass increases the higher 
Et cuts yield better expected limits. 

We use a Bayesian method Q to calculate cross sec- 
tion limits. Inputs to the calculation are signal accep- 
tance, estimated background, and observed data. The 
signal acceptance and background are assigned priors 
and modeled via a Monte Carlo method that allows cor- 
relation of uncertainties between acceptance and back- 
ground. In our analysis, the JES and luminosity un- 
certainties in the acceptance and in the background are 
correlated. The expected limits are calculated by simu- 
lating observed data based on the expected background 
with Poisson fluctuations. 

Figure [2] shows typical invariant mass distributions 
reconstructed for each particle type for a mass of 600 
GeV/c 2 using the optimal set of Et cuts in each case. 
The WW invariant mass distributions are shown for G* 
and Z', and the W Z invariant mass distribution is shown 
for W'. The optimal set of E T cuts for G* at 600 GeV/c 2 
is from the second series with Et > 120 GeV, while both 
Z' and W' favor the first series with E T > 60 GeV. The 
background compositions, as shown in Table HI are found 
to be more sensitive to the different sets of Et cuts than 
to the different decay types (WW or WZ). For Z' and 
W', the QCD jet background has a much lower contribu- 
tion owing to the stricter Et requirements. 

Without a statistically significant excess above the ex- 
pected background in the invariant mass plots, we calcu- 



TABLE I: Percentage fractional background compositions in 
Fig. [2] The uncertainties include both statistical and system- 
atic uncertainties. 





G*(WW) 


Z'(WW) 


W'(WZ) 


W+jets 


31.8 ± 8.2 


33.0 ± 10.0 


36.8 ± 9.7 


it 


19.6 ± 2.7 


35.1 ± 4.0 


37.4 ± 5.2 


WW 


10.7 ± 3.2 


15.2 ± 2.8 


13.4 ± 3.2 


QCD jets 


32.7 ± 6.5 


5.1 ± 1.0 


5.6 ± 1.1 


Others 


5.3 ± 0.9 


2.4 ± 0.9 


3.4 ± 1.0 



late the cross section limits at 95% confidence level (C.L.) 
for the observed data. Figure [3] shows the observed and 
the expected 95% C.L. cross section limits overlaid with 
theoretical cross sections. The theoretical cross sections 
for G* and Z' are calculated from pythia version 6.216, 
and a constant K- factor of 1.3 is applied to take into ac- 
count the NLO correction The theoretical cross 
section for W' is derived from a NLO calculation [23j j . 
The upper right inserts in Fig.[3]show ratios of the limits 
to the theoretical cross sections. Where the ratio is below 
one the mass region is excluded. Table ILT1 summarizes the 
mass exclusion regions from the figures. 



TABLE II: Mass exclusion region at 95% C.L. with k/M Pi = 
0.1 for G*, and £ = (M w /M v ) 2 (C = 1) for Z' and W . 





G* 


Z 1 


W 1 


Expected Exclusion (GeV/c^) 
Observed Exclusion (GeV/c 2 ) 


< 632 

< 607 


257-630 
247-544 


381-421 
285-516 



The results shown in Fig. |3] and Table [TT| for Z 1 and 

W' are based on a gauge coupling mixing factor of 
£ = C ■ (M w /M v ) 2 , with C = 1. Since signal accep- 
tance is the only quantity that changes with £ in the 
cross section limit calculation, at each mass point we re- 
evaluate signal acceptances for different £ values and cal- 
culate cross section limits as a function of £. Comparing 
the calculated and theoretical cross sections as a function 
of £, a £ exclusion region is derived at each mass point. 
These Z' and W' exclusion regions are shown in Fig. [4] 
The branching ratio of Z' or W' to fermions decreases 
as £ increases. This is opposite to the diboson decay 
modes where branching ratios increase as £ increases. 
Most Z' or W' search results @, 0] report mass limits 
along the £ = (Mw / My) 2 line and we have also done so 
for comparison. However, the diboson decay modes and 
the fermionic decay modes are sensitive to different parts 
of the gauge coupling strength phase space, so searches 
for bosonic and fermionic decays of Z' and W' are com- 
plementary to each other. The W' result shown in Fig. 0] 
is significantly improved compared to the previous result 
from CDF Run I 0. The Z' result shown is the first to 
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set an exclusion region as a function of £ and mass. 

In conclusion, we have searched for new particles de- 
caying into a pair of bosons in the electron, f,T, and 
two jets final state. In data from an integrated lumi- 
nosity of 2.9 fb _1 , no significant excess over the standard 
model prediction is observed. Cross section limits at 95% 
C.L. and mass exclusion regions have been obtained for a 
Randall-Sundrum graviton, Z' and W bosons. The W 
exclusion region in the £ — Mw' plane has been extended 
significantly compared to the previous measurement. We 
have also presented the Z' exclusion region in the £ — Mz> 
plane for the first time. We set the most stringent mass 
limits on W and Z 1 bosons. 
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FIG. 2: Invariant mass distributions with optimal set of Et cuts for 600 GeV/c 2 signals. Left to right: G*(WW), Z'(WW), 
W (WZ). G* and Z' are the same decay mode (WW) but with different optimal selections. Z' and W' have the same optimal 
selection but different decay modes (WW vs WZ). Superimposed on top of the backgrounds are signals at 600 GeV/c 2 mass 
corresponding to the theoretical cross sections. 




FIG. 3: Cross section limits at 95% C.L. Left to right: G* , Z' , W' . Inserts at upper right are cross section limits divided by 
the theoretical cross sections. Yellow and green bands show la and 2a bands on the expected limits, respectively. 




FIG. 4: Z' (left) and W' (right) exclusion regions as a function of mass and £. The £ = (Mw / Mv) 2 (i.e., C = I) lines indicate 
pythia defaults and are commonly used for mass exclusion regions. The vertical lines mark the results as shown in Table ITTl 
Also shown in the W' plot is the CDF Run I result. 



